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I.  Introduction 

Motion  of  incompressible  and  gaseous  media  is  represented  by 
the  Navier-Stoke  differential  equations.  For  two-dimensional  flow  and 
constant  properties,  these  are,  in  Cartesian  coordinates: 
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In  addition  to  these,  we  have  the  continuity  equation: 
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where  u  and  v  are  the  velocities  in  x  and  v  directions,  X  and  I 
are  the  mass  forces,  p  is  the  pressure,  p  the  density,  and  p  the 
viscosity.  To  date,  a  solution  of  the  non-linear  differential  equation  (1) 
has  not  been  found;  not  even  if,  as  a  simplification,  we  restrict  ourselves 
to  stationary  properties  which,  from  a  mathematical  point  of  view,  will 
exclude  all  cases  of  turbulent  flow.  More  responsive  to  a  mathematical 
treatment  are  stationary  flows  at  very  small  ar.d  at  large  Reynolds  numbers.* 

In  the  first  case,  the  convection  terms  (u  >  v  etc.)  in  Lq.(l)  can 
be  neglected ;  and  in  the  second  case,  the  effect  of  viscosity  can  be  neglected 
as  long  as  no  flow  separation  takes  place  (potential  flow).  This  is  true 
except  for  a  very  thin  layer  directly  on  the  surface  of  the  boundary  (boundary 
layer),  in  which  the  longitudinal  velocity  decreases  with  a  very  steep  gradient 
to  aero  on  the  boundary,  and  in  which  the  visco.ity  conditions  miy  not  be 
neglected.  The  flow  in  this  boundary  layer  can  also  ie  treated  mathematically 
as  long  as  the  layer  remains  laminar.  Such  boundary  layers  develop  on  all 


*  Reynolds  number  is  understood  to  be  Ud  ,  where  U  denotes  a  characteristic 

y 

velocity,  d  a  characteristic  length,  and  the  kinematic  *xscosity. 
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bodies  subjected  to  a  flow  ani  are  related  to  the  flow  resistance  of  the 
body.  This  originates  from  the  forces  acting  on  the  body,  which  can  be  broken 
down  into  two  components:  Forces  normal  to  the  cody  (pressure  forces);  and 
those  tangential  to  the  body  (shear  stresses),  due  to  the  viscosity  of  the 
flowing  medium.  If  no  flow  separation  from  the  body  takes  place,  the  pressure 
forces  can  tt  calculated  in  a  good  approximation  from  the  potential  flow  around 
the  body;  and  the  shear  stresses,  from  the  boundary  layer  How. 

A  body  for  which,  to  begin  with,  nothing  but  shear  stresses  can  be 
expected  is  a  flat  plate  in  longitudinal  flow.  For  this,  in  190 4,  Prandtl 
was  first  to  show  that,  even  at  high  Reynolds  numbers,  the  drag  could  be 
calculated  analytically  from  the  Nsvier-Stoke  equations.  Here,  the  Reynolds 
numbers  are  assumed  to  be  high,  yet  not  too  high,  so  that  the  flow  remains 
laminar  throughout  (Re  *  —  <  300,000  ,  x*  distance  from  the  leading  edge  of 
the  plats).  Then,the  thickness  of  the  boundary  layer  normal  to  the  plate 
(y  -  direction)  is  small  as  compared  with  its  dimension  in  the  direction  of 
the  plate  (x  -  direction),  so  that,  as  was  proyen  by  Prandtl^,  we  may 

limit  ourselves  in  Eqs.  (1)  to  the  u-components ;  and  further,  is  small 

d*"u  dD  dx 

as  >  so  that  it  may  be  neglected;  and  we  may  write  ^  *  0.  Thus,  for 

constant  properties  we  obtain  (introducing  kinematic  viscosity  i>  -  ^  ): 
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with  the  boundary  con^iti.-ns  (U  constant  free-stream  velocity  in  frontof  the 
plate) 
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By  introducing  the  non-dimensional  distance  from  the  boundary  f  and  the  stream 
function  ty, 
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Equations  (3)  and  (3a)  can  be  reduced  to  conmon  differential  equations, 
where  the  continuity  condition  (3a)  is  satisfied,  by  writing  u  - 
and  .  Thus,  we  obtain 

f  +  f f  -  0  .  (6) 

A  numerical  solution  of  this  differential  equation,  (6),  was  given  by 

Blasiua.  |cf.  Refs.  3,  u\  Today,  voluminous  literature  is  available  on 

boundary  layers,  even  for  those  at  pressure  dro,_s  and  pressure  rises  (^  4  0). 

dx 

If  there  is  a  temperature  difference  between  the  body  and  the 
flowing  medium,  a  heat  exchange  between  the  two  takes  place,  for  .hich  the 
heat  conduction  equation  holds: 


dT  dT  .  dT 
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Where  T  denotes  temperature,  X  thermal  conductivity,  and  c^  specific  heat 
at  constant  pressure.  To  begin  with,  the  heat  ;  reduced  by  friction  has  not 
been  considered  in  the  differential  equation  (7).  This  is  justified  as  long 
the  velocities  encountered  are  not  toe  nigh  and  the  temperature  differences 
are  not  too  small.  All  conditions  for  the  flow  field  are  analogously  valid 
for  this  case  alsoj  in  particular,  near  the  surface  of  the  body  a  ”  Uieraal 
boundary  layer"  develops,  in  which  the  entire  temperature  drop  between  the 
body  and  the  surrounding  medium  takes  place.*’  If  the  corresponding  flow  field 
is  already  known, then  the  solution**  of  (7)  for  the  flat  plate  at  zero  angle 

/  r  \ 

•f  attack  can  be  reduced  to  a  pure  quadrature,  as  was  shown  by  E.  Polhausen: 

„_T-£  .Km  .  .  .  fl 
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where  T0 denotes  the  temperature  at  the  wall,  and  that  at  the  outer  edge  of 
the  boundary  layer***.  Obviously,  the  boundary  conditions  y  «  0,  T  «*  T  , 
y  »oo,  and  T  «  T^  are  satisfied. 

QL  m  m 

*  Therefore  rrr  may  again  be  neglected  with  respect  to  ;  furthermore, 
ojt  oy 

we  limit  ourselves  to  stationary  temperature  fields  (~»0). 

OL 

**  With  the  simplifying  assumptions  mentioned  above. 

***  Pr  is  the  Prandtl  number  ^  ,  and  f  is  defined  by  Eq.  (5). 
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If,  in  the  flow  field,  there  are  media  of  different  kinds,  the 

concentrations  of  which  have  local  differences,  then  Eq.  (7)  is  replaced 

by  a  similarly  constructed  diffusion  equation  in  which  the  temperature  T 

is  replaced  by  the  concentration  c,  and  the  thermal  conductivity  a  -  - 

pcp 

by  the  diffusion  coefficient  k . 

In  applications  to  practical  problems  of  heat  transfer  and 
diffusion,  difficulties  arise  due  to  the  fact  that  the  stipulation  of  con¬ 
stant  properties  is  frequently  not  satisfied,  since  most  properties  are 
functions  of  temperature  and  then  flow  -  and  temperature  fields  are  no 
longer  independent  of  each  other.  A  similar  difficulty  arises  in  the 
domain  of  heat  transmission.  On  the  basis  of  analogy  considerations,  we 
can  represent,  for  constant  properties,  the  results  of  heat  transfer  measure¬ 
ments  in  the  form: 


Nu 
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where  Nu  is  the  non-dimensional  Nusselt  number^ 
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\  >  T  >  T  >  etc*>  3X9  the 

proportions  determining  the  geometrical  analogy  conditions,  and  d  is  a 
characteristic  dimension  of  length.  The  question  arises'  whether  the  vari¬ 
ability  of  the  properties  could  net  be  taken  into  account  by  substituting 
into  (9)  the  properties  for  a  suitably  selected  temperature,  so  that,  again, 
all  measurements  could  be  represented  in  a  standard  fora. 


In  diffusion  problems,  the  density  may  depend  upon  the  concentra¬ 
tion  if  the  problem  involves  high  concentrations  and  media  of  different 
densities.  This  corresponds  to  variable  properties  in  temperature  problems. 
However,  another  phenomenon  ray  appear:  At  high  concentrations,  the  normal 
velocity  v(0)  no  longer  vanishes  at  the  boundary,  a  fact  pointed  out  previously 
by  Nusselt. ^  If  a  liquid  covers  a  wall  and  is  made  to  evaporate  by  a  gas 
flowing  past  this  wall,  then  some  liquid  enters  continously  into  the  flow 
field  and  we  obtain,  therefore,  v(0)  >  0.  If  vapor  condenses  on  a  wall,  or 
if  a  gas  is  bound  to  it  by  a  chemical  reaction,  then  v(0)  <0. 


( 
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Presently  there  are  two  papers  available  on  the  laminar  boundary 
layer  on  a  flat  plate  in  parallel  flow  with  variable  properties  -  one  by 
L.  Crocco^*^  and  the  other  by  V.  Karman  and  Tsien.^  In  both  papers,  the 
differential  equations  are  changed,  by  means  of  a  transformation  to  new  vari¬ 
ables,  into  a  form  that  is  different  from  the  common  boundary  layer-  equations. 
Crocco  obtains  two  simultaneous,  second  degree,  differential  equations  which 
he  solves  for  a  gas  with  the  Prandtl  number  0.725  (Air),  Karman  and  T3ien 
examine  the  special  case  of  Pr  »  1  and  have  to  solve  but  one  differential 
equation,  since,  in  this  case,  the  temperature  field  has  a  simple  relation 
to  the  velocity  field. 

This  study  begins  with  an  application  of  a  method  of  solution  (to 
which  little  reference  is  made  in  the  literature)  for  the  flow  with  constant 
properties  over  a  flat  plate.  This  method  of  solution  was  suggested  origi¬ 
nally  by  Piercy  and  Preston,  employing  Polhausen’s  formula  (8).  Futhermore, 
it  is  shown  how  the  boundary  layer  flow  with  variable  properties  over  the 
flat  plate  can  be  calculated  simply  and  clearly  by  means  of  this  solution 
method;  and  also,  in  diffusion  problems,  how  the  finite  normal  velocity  at 
the  boundary  can  be  taken  into  account. 

Let  us  limit  ourselves  to  laminar  boundary  layers,  since  analytic 

* 

treatment  of  the  corresponding  turbulent  case  is  substantially  more  difficult 
and  complicated.  Laminar  boundary  layers  appear  on  bodies  in  the  vicinity  of 
the  stagnation  point  and  extend  from  there  over  a  distance  that  increases 
with  decreasing  Reynolds  number,  it  being  stipulated,  of  course,  that  no 
flow  separation  occurs.  Such  conditions  present  themselves  particularly 
(1)  in  the  flow  of  viscous  fluids;  (2)  in  cases  in  wnich  the  body  in  the  flow 
is  very  small  (e.g.,  particles  in  coal  dust  firing);  and  (3)  in  diffusion 
problems  where,  frequently,  only  low  velocities  appear. 

Primarily,  this  study  aims  at  fundamental  information,  while 
practical  applications  are  possibly  to  be  treated  in  a  later  report. 
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II.  Solution  of  the  Boundary  Layer  Equation  for  Variable  Properties 


Taking  variable  properties  into  account,  the  boundary  layer 
equations  for  the  velocity  and  temperature  fields  on  the  flat  plate 
arJ5> 
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The  notations  are  the  same  as  in  (1),  (2),  and  (7).  Again ^ frictional 
heat  is  as  yet  neglected. 

For  constant  properties,  (10)  and  (11)  can  be  reduced  to  a 
common  differential  equation  (cf.  introduction  and  Eqs.  (3)  through  (6)  ) 
by  the  assumption  that  u  and  T  are  a  function  of  but  one  (non-dimen¬ 
sional)  coordinate,  5  “  2  'fix  *  S*nc®  t^le  Properties  depend  only  upon 
temperature,  the  conclusion  is  obvious  that  the  same  simplification  is 
possible  for  variable  properties  also.  Wo  write 


u 

tr 
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where  U  is  the  velocity  at  the  outside  of  the  boundary  layer,  and  Tfi  and 
T-^  are  the  wall  temperature  anJ  th<.  temperature  at  the  outer  limit  of  the 
boundary  layer,  respectively.  The  magnitude  2/^  in  the  non-dimensional  $ 
denotes  the  kinematic  viscosity  at  the  fixed  temperature  1^,  for  which 
the  wall  temperature  (k  *  0),  or  the  temperature  at  the  outer  limit  of  the 
boundary  layer  (k  *  1),  is  a  suitable  selection.  The  boundary  conditions 
for  the  flow  and  temperature  fie 

y  -  0,  l  -  0  ,  oo"  0  , 

y»oo,  l  m  00  ,  (J  -  1  , 
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Furthermore,  let 

£  -  °  -  V(5)  ,  r  -  x(5)  ,  (H) 

Hi  K 

where  the  index  k  denotes  the  'property  at  the  temperature  . 

The  specific  heat  Cp  is  assumed  to  be  constant,  which  is  true 
for  most  gases  and  fluids.  To  begin  with,  we  obtain  from  Eq.  (10a) 

pV  “Vttx  (PUS~  fpud$)>  (15) 
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which  leads  from  (10)  and,  after  (14)  and  (15)  have  been  introduced,  to  a 
fora  that  is  suitable  for  further  calculations: 
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From  (16),  the  following  expression  can  be  derived  for  cj  ,  if,  temporarily, 
(16)  is  considered  as  a  differential  equation  for  /  and,  in  addition, 
if  f  is  regarded  as  a  known  function  of  3  : 

.  -/  k  ** 
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where  the  constant  of  integration  is  given  by  the  boundary  condition  (13). 
Similarly,  we  obtain  for  the  non-dimensional  temperature  8  the  expression 


-Pr, 
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(18) 


where  Pr, 


a  ^  cpk 

°k  \ 

at  the  temperature  T^, 


denotes  the  Prandtl  number  with  the  prope-^ies 


—8— 


For  constant  properties  the  velocity  and  tempera¬ 

ture  fields  are  independent  of  each  other,  and  we  obtain  from  (18)  the  solu¬ 
tion  given  by  Polhausen^  for  the  temperature  field  (cf.  Eq.  8  )  which 
represents  an  integral  equation  for  the  velocity  field  for  Pr  *  1  -  a), 

as  may  be  found  by  a  comparison  of  (17)  and  (18).  While,  for  a  known  velocity 
field,  solution  for  the  temperature  field  according  to  (8)  or  (18)  is  possible 
by  means  of  simple  quadrature,  the  Calculation  of  the  velocity  field  involves 
the  difficulty  that,  in  (17) ,  the  still  unknown  velocity  appears  on  the  right- 
hand  side  in  the  expression  for  f.  In  the  method  of  solution  given  by  Piercy 
and  Preston,  we  proceed  from  an  arbitrary  approximation  for  a)  which  is  then 
employed  to  calculate  f  according  to  (16),  upon  which  J(|)  is  calculated  by 
means  of  (17).  In  this  process  it  should  be  borne  in  mind  that  constant  proper¬ 
ties  are  stipulated,  and^thus,  that  /  *  o  1.  In  this  manner  an  improved 

value  of  cj  is  obtained  which  represents  the  initial  value  for  a  new  calculation, 
etc.  Figure  1  shows  graphically  the  individual  steps  of  the  approximation.  As 
initial  solution  c/°  ,  the  intentionally  inaccurate  approximation  cJ  «  1  over 

the  entire  boundary  layer  has  been  selected;  the  associated  first  approximation 
is  given  by  the  error  integral.  After  the  third  approximation,  the  shear 
stress  on  the  wall  deviated  but  4.5  percent  from  the- exact  value.  Instead  of 
continuing  this  process  mechanically,  the  expected  final  solution  was  estimated 
from  the  course  of  the  preceding  approximations  and  then  used  as  a  basis  for 
the  subsequent  step  in  the  approximation.  In  this  manner  the  solution  o  was 
obtained  with  only  l/2-percent  error  in  the  shear  stress. 

For  this  method  of  solution,  the  improvement  obtained  by  oach  approxi¬ 
mation  step  can  be  estimated  quantitatively:  Eqs.  (17)  and  (18)  are  identical 
for  constant  properties  and  Pr  ■  1.  Let  us  assume  that  we  dealt  with  an 
approximate  solution  of  3uch  a  nature  that  the  associated  I  coordinate 
differed  by  a  constant  factor  Y  from  the  5  coordinate  of  the  exact  solutions. 
Then,  obviously,  we  have  also  (f)^  -  y  f,  a  comparison  of  (17)  and  (18) 
yields  that  the  influence  of  the  factor  V  is  equal  to  that  of  the  quantity  Pr 
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on  the  temperature  field.  Polhausen^  found,  or.  tne  basis  of  his  numerical 

^  - 

calculations,  that  the  heat  transfer  coefficient  is  proportional  to  ; 

consequently,  using  the  netnod  mentioned  abGve,  the  shear  stress  at  the  wall 
is  subject  to  an  error  in  every  approximation  step  that  is  only  about  one-third 
of  the  error  of  the  preceding  one. 

In  the  case  of  variable  properties,  these  solution  steps  of  a 
11  mathematical"  nature  can  be  combined  with  steps  of  a  "  physical"  nature: 

Step  1 

As  the  initial  point,  the  known  solutions  for  constant 
properties  are  assumed. 

(a)  For  the  velocity  profile,  the  solution  by  Blasius. 

(b)  For  the  temperature  field,  the  solution  by  S.  Polhausen. 

Step  2 

(a)  Calculation  of  the  velocity  profile  according  to  (17), in 
which  the  temperature  profile  from  3tep  1  (b)  accounts 
for  the  temperature  dependence  of  the  properties. 

(b)  Calculation  of  the  temperature  field  according  to  (18) 

by  means  of  the  velocity  profile  of  step  2  (a);  dependence 
of  tfte  properties  upon  temperature  is  taxen  into  account 
in  the  same  manner  as  in  step  2  (a). 

This  procedure  is  repeated  until  the  final  solution  is  sufficiently  exact. 
Generally,  3  to  4  repetitions  will  yield  satisfactory  results. 

First,  let  us  elaborate  on  tr.e  effect  of  the  dependence  of  the 
properties  on  temperature.  Of  special  interest  in  the  velocity  or  tempera¬ 
ture  field  are,  chiefly,  single  values,  suer,  as  the  velocity  gradient  at  the 
vail  (for  the  calculation  of  the  shear  stress;  or  t.ne  heat  transfer  coefficient. 
The  idea  suggests  itself  to  take  into  account  tr.e  temperature  dependence  of  the 
properties  by  substituting  the  properties  at  a  suitably  selected  ten.  erature 
into  the  ** isothermal11  formulas  (i.e.,  the  formulas  for  constant  properties). 

If  ve  select  as  reference  temperature  that  at  one  o:  the  two  limits  of  the 
boundary  layer,  then,  on  the  basis  of  physical  aspects  as  *ell  ac  on  the  basis 
Of  the  formulas,  it  will  be  found  that,  within  the  boundary  layer,  an  increase 
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in  viscosity  or  density  as  compared  with  the  values  at  the  limits  of  the 
boundary  layer  will  result  in  a  drag  increase  a c  compared  with  the  iso- 
thernal  flow.  Similarly,  an  increase  ir.  tnermal  conductivity  and  density 
results  in  increased  thermal  dissipation.  However,  the  question.  How 
great  is  the  effect  of  the  unconstancy  of  the  properties?  depends  upon  the 
ratio  of  the  boundary  layer  thickness,  of  the  temperature ,  ana  of  the  velc 
city  field.* 

Let  us  illustrate  this  by  means  of  the  following  case,  which  has 
also  some  practical  significance:  The  thermal  boundary  layer  is  assumed  to 
be  very  small  *s  compared  ..ith  the  flow  boundary  layer,  a  condition  encountered 
at  large  Prandtl  numbers  (viscous  fluids).  In  this  case,  then,  the  variation 
of  t;.e  properties  within  the  thermal  boundary  layer  obviously  may  be  neglected 
for  the  determination  of  the  shear  stress,  which  will  be  found  to  be  the  same 
as  if  the  temperature  at  t  he  outer  edge  of  the  boundary  layer  would  extend  to 
the  wall.  The  same  holds  also  for  the  velocity  profile,  with  the  exception 
of  tne  sma-i  region  w_i.tr.in  the  thermal  boundary  layer  in  which  the  velocity 
profile  is  deformed  in  accordance  with  the  viscosity  variation.  For  the 
temperature  profile,  however,  exactly  this  region  is  decisive.  For  the 
velocity  gradient  on  the  wall,  we  obtain,  from  the  equality  of  the  shear 
stress. 


(19) 


where  t..e  Index  . »  denotes  t..e  "  isothermal”  flow  with  the  properties  at 
the  temperature  T^.  In  trJLs  case,  variability  of  density  has  no  effect  on 
the  flow  field.  This  can  also  be  derived  mathematically  from  Eqs.  (17)  and 
^18).  The  conditions  for  the  temperature  field  are  discussed  in  the  following 
paragraph  by  means  of  two  examples. 


*  According  to  E.  Polhausen,  the  ratio  of  the  two  is  proportional  to 

ViT 


IH.  Flew  and  Temperature  Fields  in  Viscous  Fi uids 


In  accordance  vitr.  the  nature  of  the  properties  of  viscous  fluids, 

the  velocity  and  temperature  fields  were  calculated  with  the  assumption  that 

only  viscosity  varies  with  temperature  according  to  the  following  formula: 

{  T,  ♦  T  b 
u  /  x  c 

M*'(j  *TC  1 

where  b  and  T  are  constants  vhicn  are  to  be  selected  in  such  a  manner 
c 

that  the  temperature  dependence  is  as  well  represented  as  pc  £  si  hie.  Let 

the  index  k  have  the  value  0  or  1,  depending  cn  the  selection  of  the  non- 

dimensional  5  .  Here,  b  «  3  was  selected  (e.g.,  vise: us  1-tricutir.g  oil,, 

r-Q  p 

and  the  two  cases  of  a  heated  and  coded  plate  were  culculatec  with  — *  *  — 

o 

and  S,  respectively,  and  PrQ  -11.5  and  IX  .  Thus,  tre  two  exanples^deal 
with  the  sane  fluid  and  involve  equal  temperature  differences,  £ ir.ee,  in. 
both  cases,  the  respective  Ft  is  forced  -  itr.  t  e  properties  of  tne  wail  temp— 
e nature.  If  I  is  selected  as  reference  temperature ,  we  octain  frem  .X, 


if-? 


-i  ■  - 1  ! 

V/  ^  /  / 


The  result  of  the  calculations  accord  u  to  the  iterative  method  o. 
the  preceding  paragraph  nay  ce  recognized  in  rigs.  ana  }.  In  both  graphs 
the  non-dimensional  distances  from  the  wall,  and  5^  ,  forced  with  and  ^  , 
are  plotted  on  the  abscissa  at  the  scale  1:  (Fig.  2)  and  -*8  :  1  '(Fig.  ; , , 

so  that  but  one  point  on  the  acscissa  correspond*  to  every  oust once  from  the 
wall  y  whether  this  is  determined  fc%  the  notation  ^  or  .  Thus.  ir. 
both  graphs  the  location  ox  the  curves  vitr.  re-pect  to  eacr.  other  :c. — *  ;  .:,.s 
to  the  real  conditions.  Besides  the  solution  u,  for  th»*  -  ,icr. 

three  steps  were  required,  the  isothermal  velocity  p-’-'l^es  (—70  a^nd  for 

constant  properties  at  the  temperatures  and  ,  respectively ,  were  plotted 

v  >  P  1 

as  a  function  of  the  non-dimensional  .oordinates  *  and  5.  .  For  the  isotneruai 

o  1 

teeperature  profiles  (0)o  and  (8)1#  the  Pr  numbers  also  «ere  substituted  at 


— vi-c— 


tee  temperatures  ?  and  T, .  For  example,  in  Figure  2  ?r  -  12,5  *nd 

o  1  o 

Pr,  “  IX.  In  the  following  table,  'J'  and  =<  denote__the  shear  stress 


at  the  vail  ana  neat  transfer  coefficient  j*<-  '*  ^To4  aXiJ* 

are  t..e  corresponding  values  at  isot hemal  flo-  with  the  viscosity  at  the 
temperature  T,  and  T  ,  respectively.  The  sane  holds  for  (TT),  and  . 

♦  Q  O  y  X 

Table  1 


^o 

?r 

Vo 

(r0)x 

Heated  wall 

0.125 

12.5 

0.841 

2.38 

Cooled  wall 

8,0 00 

100.0 

1.05 

0.382 

/dcp\  /de 

Woi  (d$e 


Although  in  these  case-  the  thickness  of  the  thermal  boundary  layer  is  by  no 
means  negligible  as  compared  with  the  flow  boundary  layer,  the  shear  stress 
can  be  calculated  very  well  with  the  aid  of  the  isothermal  formula  with  the 
viscosity  at  the  temperature  ?^.*  The  stipulations  of  Sc.  (19)  for  Pr  >  10 
are  thus  given. 

The  conditions  for  the  heat  transfer  coefficient  are  more 
involved.  From  (S),  it  follows  that  the  heat  trarr  fer  coefficient  is 
proportional  to  g(? r)  ,  w'-ere,  according  to  S.  Polhausen,  g  ■  0.664  ^/?r 

is  true  with  r.irh  accuracy.  Considering  that  Pr  -  ^  ,  it  follows  that  the 
heat  transfer  coefficient  i3  inversely  proportional  to  the  sixth  root  of  the 
viscosity.  7nus,  the  ratio  of  heat  transfer  coefficients  fa)Q  to  (=C)^  is 

Wo  »/^lV  (22) 

J&1  \»0j  * 

r  second  inuicati^n  for  tr.e  neat  transfer  coefficient  is  oct-iined  from  the 
velocity  gradient  a*  the  waul.  From  Eq.  (19),  with  consideration  of  (12), 
it  follows  (if  u.e  indices  at  the  brackets,  i^  and  ,  denote  isothermal 

flow  at  the  temperature  Tq  and  T^  ,  respectivexy) : 


*  Refer  to  column 


in  Table  1  ani  cohere  it  with  column 


iumn 


o 
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These  relations  also  are  c-nfimed  by  Figures  2  and  3- 


(23) 


Obviously,  by  c>>ans  of  these  fonsuias,  limits  c an  be  given  for  the 
heat  transfer  coefficient  (also  compare  Figs.  2  and  3),  which  we  shall 
express  in  terns  of  multiples  of  the  heat  transfer  coefficient  .  une 


6  ' 


liait  is  given  by  (<)^  ■  /'  (°<)o  (according  to  £q.  ^2),  because  the 

^  ^1 

velocity  profile  (cj)^ .  yields  higher  velocities  on  all  points  with  a  cooled 
plate  and  lower  velocities  with  a  heated  plate.  The  other  linii  is  given  by 
a  velocity  profile  of  isothermal  fora  for  vrJ.cn  the  abscissa  sc«^e  is  changed 
by  a  factor  y  in  such  a  manner  that  its  gradient  at  the  wall  coincides  with 
that  of  the  actual  velocity  distribution.  Since  the  limit  mentioned  above 
for  the  heat  transfer  coefficient  is  to  be  expressed  by  (x)q  nultiplied  by 
a  factor,  y  is  to  be  calculated  from  the  ratio  of  the  true  velocity  gradient 
to  that  at  the  temperature  Tq  (second  part  of  Eq.  a3).  From  the  note:  on  the 
convergence  of  the  solution  method  in  paragraph  II,  it  then  follows  tr.^t  the 
heat  transfer  coefficient  is  proportional  to  -y  y  and  we  obtain  as  the  limit 

-\  fi*. 1  (ol)  .  Summarizing,  we  obtain  for  the  limits  ol  tne  heat  transfer 

V  ]T 

coefficient,  if  only  viscosity  varies. 


(24) 


where  the  upper  signs  stand  for  t  he  case  of  the  heated  plate  and  ti.e  lo.-.er  signs, 
f®r  the  cooled  plate.  Accordingly,  in  agreement  with  the  special  examples  in 
Table  1,  the  heat  transfer  coefficient  (<*)q  approaches  closest  the  true  heat 
transfer  coefficient  oC. 


Thus,  for  viscous  fluids  (Pr  >  10)  for  which  variation  of  the 
viscosity  is  of  primary  cons  iteration,  the  following  approximate  rule  is 
obtained:  In  the  calculation  of  the  drag,  the  properties  are  to  be  refer¬ 
red  to  the  temperature  at  the  outer  Unit  of  the  boundary  layer;  in  the 
calculation  of  the  heat  transfer,  the  properties  are  to  be  referred  to  the 
wall  temperature. 


IV.  The  Flow  and  Temperature  Fields  for  tne  Case  Pr  ■  0.7(Alr) ,  r.'nere  All 
Properties  Are  Functions  of  Tenperature 

r*  o  x 

In  tr.e  tenperature  range  -50  to  140  C,  the  properties  of  air  can 
be  represented  by  the  following  formulas: 


..  m  0.7S0  ■  -  -1  . 

p  T  ;p"K2t  ;  x 

where  T  is  tenperature  in  absolute  degrees. 

Yt'iih  T  -  T 

o  1 


*3 


.0.821 


T1 


we  obtain 


p  0.78 

-  -  i  -  1  ♦  0(1  -  8) 


and  similar  expressions  for  o  and  ■%. 

The  calculation  according  to  the  method  given  in  paragraph  II  was 
carried  out  for  a  heated  plate  with  the  values  of  0  *  1/4  and  1/2  and,  in  the 
velocity  and  tenperature  fields,  yielded  only  moderate  deviations  from  the 
form  at  isothermal  flow  (Table  2).  In  order  to  calculate  the  conditions  at 
greater  temperature  differences,  the  case  of  *  20°C  and  Tq  “  620°C  (0  **  2.05) 
was  calculated.  According  to  Figure  4,  the  velocity  and  temperature  fields 
show  a  substantial  variation  of  the  fora  at  constant  properties.  Again, 
and  ^  are  formed  with  tne  properties  at  tne  temperatures  Tq  and  T^,  respec¬ 
tively.  For  both  fields,  a  substantial  increase  in  trie  boundary  layer  thick¬ 
ness  is  obtained.  Nevertheless,  the  values  of  the  wall  shear  stress  and  of 
the  heat  transfer  coefficient  show  but  little  deviation  from  the  values  at  con¬ 
stant  properties,  where  it  is  unimportant  whether  they  are  referred  to  the  temp¬ 
erature  at  the  wall  or  to  that  at  the  outer  limit  of  the  boundary  layer.  In  the 
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case  of  air,  this  car  be  explained  by  tne  fact  tnat^wit  in  tne  boundary 
layer,  the  increase  in  viscosity  with  increasing  temperature  acts  as  a 
drag  increase,  tne  decrease  in  density  act?  as  a  drag  decrease,  and  both 
effects  practically  eliminate  each  other  at  Pr  »  0.7,  where  thermal  and 
flow  boundary  layers  are  approximately  equal  in  size..  The  conditions  for 
the  temperature  field  are  almost  exactly  alike  because  thermal  conductivity 
is  related  with  temperature,  as  is  viscosity. 

Table  2 


{ 


3% 


c 


i  * 


c- 


Keating 

fa 

<0 

fa 

oL 

(fa 

'  Jo 

_ o 

H 

CDfjun 

i'—'i 

w. 

0  •  1/4 

1.00 

1.02 

1.00 

1.01 

0.436 

0.556 

0  -  1/2 

1.00 

1.05 

1.00 

1.02 

0.420 

0.485 

T  -  620° C 

0 

Tx«  20° C 

0.93 

1.11 

0.96 

1.03 

0.235 

0.286 

In  similar  manner,  frictional  heating  coin  be  taken  into  account, 
in  wuich  case  a  tern  p  (dy)  aus^  be  added  on  the  right  side  of  Eq.  (11), 


and  the  solution  becomes 

6-iARM 

^ 


-  I  ^  e 


dS 


„  Prk  &  Te  Pi  i  (5  .  .  ,.*■  -  R(«)  -  m)  . 

B(5)  -  2  \  x  1  J  ?<*>  e  '  dS!  e  VJ/  d<  t  (25) 


H($>  -  Prw  <»5  , 


L 


«.  -  £ 
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Here,  also,  the  iteration  method  can  be  carried  out  even  though  the 
calculation  work  involved  is  somewhat  greater.  For  constant  properties, 
Eq.  (25)  reduces  to  the  solution  given  by  E.  Eckert. By  a  suitable 
selection  of  the  constant  of  integration  in  (25),  even  the  thermometer 
proolea  (vanishing  temperature  gradient  at  the  wall)  can  be  solved. 

Considering  Crocco’s  previous  calculations  for  a  gas  with 
Pr  ■  0.725,  a  numerical  example  worked  out  by  this  new  method  has  been 
omitted  here. 


V.  Application  to  a  Diffusion  Problem 


dc  dc  ,  d  c 
u  —  ♦  v  —  -  k  t— 7 
dx  dy  dy1 


The  concentration  field  in  the  diffusion  problem  for  the  flat 
plate  can  be  treated  in  the  same  manner  as  the  temperature  field.  The 
differential  equation  is 

A*- 

>  .  (26) 

where  k  is  the  diffusion  coefficient  and  c  is  the  concentration  that 
is  defined  as  a  quantity  of  gas  or  vapor  in  unit  volume.  Let  us,  in  tnis 
case,  consider  the  properties  to  be  constant,  and,  in  particular,  let  us 
assume  that  the  densities  of  the  two  interdiffusing  media  are  approximately 
equal  in  the  entire  flow  field.  On  the  other  hand,  let  us  take  into  account 
that,  at  higher  concentrations,  the  velocity  v  at  the  wall  no  longer  van¬ 
ishes,  a  fact  pointed  out  previously  by  Nusselt.^  When  fluid  evaporates 
from  a  wall,  say,  because  of  a  gas  (e.g.,  air)  flowing  past  a  wetted  wall, 
then  some  substance  enters  continuously  into  the  flow;  therefore,  we  obtain 
at  the  wall  v(0)  >  0.  If,  on  the  other  hand,  the  vapor  from  a  vapor  gas 
mixture  condenses  on  the  wall,  or  if,  for  example,  ammonia  containing  air 
flows  past  a  blotting  paper  saturated  with  hydrochloric  acid,  then  we 


*  E.  Eckert  reported  on  a  solution  of  this  problem  at  a  meeting  of  the  VDI- 
Ausschus3  fur  Warmeforschung,  1943,  in  Bayreuth, where  an  a oproximation 
method  was  applied  similar  to  K.  Polhausen's  method  for  the  flow  boundary 
layer . 

t  An  estimate  for  tne  problem  at  hand  has  been  given  by  G.  Damkohler, 

Z.  Elektrochem,  178  (1942). 
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obtain  v(0)  <  0.  According  to  Eq«.  (100)  and  (101)  by  liusselt,^  the 
boundary  condition  for  v  is 


k  ft c ’•  1 

'  — 7  ' v(0) 

1  n 


where  c  denotes  the  concentration  of  the  gas  or  vapor  for  which  the  wall  is 
permeable;  cq,  the  corresponding  concentration  at  the  wall;  pQ  ,  the  asso¬ 
ciated  partial  pressure;  and  p, the  total  pressure. 

For  the  velocity  profile,  the  solution  of  Biasius^  is  now  no 
longer  obtained,  because  v(G)  /  0  (see  Fig.  1),  but  is  a  family  of  profiles, 
depending  on  the  value  of  v(C).  By  introducing  the  free-strean  velocity  U 
and  the  non-dimensional  5  >  we  obtain 


/ v  I  j>_  1  k  C1  ~  co  /  dC  'j  —  /[}  f 

co  /£  _  L  id5y0  v>'x 

lpo  1  c 


0  -  c 


ci  “  co 


where  c  and  c,  denote  the  concentrations  at  the  wall  and  at  the  outer 
o  1 

limit  of  the  boundary  layer.  Hence,  similar  to  (15), 


H  -  -  - 


c  -  c 
1  o 


"  ‘"(I.-1) 


First,  let  us  state  the  solution  for  the  concentration  field  wr.ich  is, 
analogously  to  (18), 

nil  v  d  -  i  f < f<5>  -«)  d5  „ 


f(5)  -  *  j  §  di . 


where  ^  is  a  value  analogous  to  the  Pr  number.  In  order  to  obtain  the 

velocity  CO  -  ,  we  simply  write  in  Eq.  (33),  *r  ■  1.  Calculation  of 

U  K 
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the  velocity  and  concentration  fields  is  acconplished  by  an  iteration 
method  similar  to  that  in  paragraph  II. 

K  incorporates  the  concentration  gradient  at  the  wall;  how¬ 
ever,  (30)  can  be  solved  for  any  K  values  and  we  can  then  evaluate 


N 


k 


c.  -  c 
1  o 

'•(I.  -  ') 


by  means  of  the  value  for 


(si 


obtained  from  the 


solution.  The  velocity  and  concentration  fields  for  the  evaluated  M  and 
N  values  nay  be  recognized  in  Figs.  5  and  6,  and  the  concentration  gradient 
at  the  wall  in  Fig.  7.  K  >  0  denotes  evaporation  from  the  wall;  M  <  Oj 
condensation  or  absorption  on  the  plate.  For  the  ratio  £  ,  the  value  0.6 
was  selected,  which  holds,  as  a  good  approximation,  for  the  diffusion  of 
water  vapor  and  aniaonia  in  air.  From  a  rigorous  point  of  view,  at  the 

higher  concentrations  stipulated,  density  and  viscosity  of  the  mixture  of 
the  two  substances  are  dependent  upon  the  concentration,  and  the  diffusion 
coefficient,  upon  the  temperature.  With  the  aid  of  the  method  described, 
such  cases  can  also  be  calculated.  If,  besides  the  diffusion,  some  heat 
transfer  takes  place,  then  the  solution  for  the  concentration  field  can 
also  be  applied  to  the  temperature  field^as  a  good  approximation.  Similarly, 
the  solution  for  the  concentration  field  presents  information  on  the  heat 
transfer  if  air  at  the  plate  is  blown  out  or  sucked  off  through,  for  example, 
a  porous  wall,  with  normal  velocities  at  the  wall  corresponding  with  Eq.  (28). 


VI.  Summary 

Following  closely  E.  Polhausen's  solution  for  the  laminar  temper¬ 
ature  field  at  the  flat  plate  in  longitudinal  flow,  formulas  are  derived 
which  permit  calculation  of  the  velocity  and  temperature  field  for  variable 
properties  by  means  of  an  integral  equation  and  an  iteration  method  based 
on  this  equation.  Accordingly,  the  following  cases  were  solved:  By  assuming 
that  only  viscosity  varies  with  temperature  and  that  the  remaining  properties 
are  constant,  the  velocity  and  temperature  fields  were  calculated  for  the 
Pr  numbeis L<.5  and  100  (viscous  fluids)  at  heated  and  cooled  plate  conditions. 
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A  more  rigorous  investigation  of  these  two  cases  yielded  an  expression  of 
general  validity.  For  a  gas  with  the  Pr  number  0.7(air),  calculations  were 
carried  out  based  on  the  assumption  that  all  properties  vary  *ith  temperature 
and  that  velocities  are  not  too  great,  so  that  frictional  heat  may  be  neg¬ 
lected.  An  increase  of  the  thickness  of  the  boundary  layers  was  obtained 
without  a  substantial  change,  however,  of  the  shear  stress  or  the  heat  trans¬ 
fer  coefficient  as  compared  uith  those  values  that  were  calculated  with  1  he 
formulas  for  constant  properties.  In  the  course  of  this  investigation,  it  was 
found  that  the  influences  of  density  and  viscosity ^and  of  density  and  thermal 
conductivity  in  the  velocity  and  temperatm e  field^on  the  wall  shear  stress 
and  heat  transfer  coefficient  are  opposed  t v  eac  .  other  and  that  they  practi¬ 
cally  eliminate  one  another.  Formulas  which  also  take  into  account  frictional 
heating  were  given,  but  elaboration  of  the  associated  calculations  was  omitted 
in  consideration  of  the  previous  results  of  the  studies  by  Crocco.  Finally, 
the  solution  methods  developed  here  were  applied  also  to  the  case  of  diffusion 
of  additional  substances,  where,  at  higher  concentrations,  finite  normal 
velocities  appear  at  the  wall,  wt.icn  results  in  a  substantial  change  of  the 
velocity  and  temperature  field. 


GRAPHICAL  REPRESENTATION  Of  THE  INDIVIDUAL  APPROXIMATIONS  IN  THE  CALCULATION 
Of  THE  VELOCITY  DISTRIBUTION  ON  THE  fLAT  PLATE  ACCORDING  TO  THE 
METHOD  PRESTON  AND  P1ERCY  (CONSTANT  PROPERTIES) 

FIG.  I 


VELOCITY  and  temperature  distribution  on  a  heateo  plate  at  variable  viscosity. 

VISCOSITY  EXPONENT  b-3.  (0),,  ANO  (<*»),,  (0),  ARE  THE  ‘ISOTHERMAL*  VELOCITY 

AND  TEMPERATURE  DISTRIBUTOR;  AND  *,  DENOTE  KINEMATIC  VISCOSITY 
AT  THE  WALL  TEMPERATURE  T0  ANO  AT  THE  TEMPERATURE  T,  AT 
THE  OUTER  LIMIT  OF  THE  BOUNDARY  LAYER 

FIG,  £ 


VELOCITY  ANO  TEMPERATURE  DISTRIBUTIONS  ON  THE  COOLED  PLATE 
NOTATIONS  ARE  THE  SAME  AS  IN  FIG.  2 

FIG.  3 


VELOCITY  AND  TEMPERATURE  DISTRIBUTION  ON  A  HEATEO  PLATE  FOR 
Pm0.7 (AIR);  ALL  PROPERTIES  ARE  VARIABLE  WITH  TEMPERATURE 

FIG.  4 
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VELOCITY  FIELD  FOR  DIFFUSION  AT  HIGHER  CONCENTRATIONS  IN  WHICH  FINITE  NORMAL. 
VELOCITIES  APPEAR  ON  THE  WALL  (COMPARE  TEXT  TO  EQS.  2S  TO  28) 

FIG.  5 


0  12  3  5 


CONCENTRATION  DISTRIBUTION  FOR  FIG.  5 
FIG.  6 


CONCENTRATION  GRADIENT  AT  THE  WALL  AND  THE  VALUE  M  AS  A 
FUNCTION  OF  N  (SUPPLEMENT  TO  Fig.  6) 


FIG.  7 
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